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Abstract. - In the spinless fermion model on a square lattice with infinite nearest-neighbor 
repulsion, holes doped into the half-filled ordered state form stripes which, at low doping, are 
stable against phase separation into an ordered state and a hole-rich metal. Here we consider 
transport of additional holes along these stripes. The motion of a single hole on a stripe is 
mapped to a one- dimensional problem, a variational wavefunction is constructed and the energy 
spectrum is calculated and compared to energies obtained by exact diagonalization. 



Introduction. - A spinless square lattice fermion model with nearest-neighbor repulsion, 
we suggest, may be a fruitful and comparatively tractable analog of the Hubbard model that 
retains many of the its properties. In particular, in the limit of infinite repulsion, near half- 
filling, the equilibrium state appears to be an array of charged, antiphase domain wall stripes 
that are stabilized by kinetic energy QJI2]- This is reminiscent of a state of charged stripes 
that has been observed in cuprates [3j and has been discussed in the spinfull Hubbard and t-J 
models 0]. Stripes in the spinless fermion model have been studied by exact diagonalization 
(ED) in ref. 2_. Here we are interested in finding to what extent a hole, moving along a 
stripe, is decoupled from the state of the stripe, and whether or not the hole tends to bind to 
a "kink" on the stripe. To this end we extend an exact mapping of a stripe to one dimension, 
introduced in refs. . Using this mapping, we construct a variational wavefunction and 
calculate the energy spectrum. 

We do not know of an existing experimental system that realizes our model, but recent 
progress in the manipulation of ultracold bosonic atoms on an optical lattice 5 , and in cooling 
fermionic atoms below degeneracy temperature (HI , leads us to expect that such a model may 
be realized experimentally in the near future. A corresponding Hard-core boson model with 
large nearest neighbor repulsion can be realized, on a triangular lattice, in adsorption of 4 Hc 
to graphite sheets 0. The phase diagram of a boson model with finite nearest neighbor and 
next nearest neighbor repulsion was studied in ref. [S]. 
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Model. - We consider spinless fermions on a square lattice with Hamiltonian 

(»i> W) 

where cj and Cj are creation and annihilation operators at site i, respectively, and (ij) means 
nearest neighbors. In this paper, we consider only the limit of zero temperature and V/\t\ = oo, 
so that neighboring fermions are forbidden, and t = 1 is the only energy scale. The maximum 
allowed filling fraction is n = 1/2, where there are two possible "checkerboard" states 9 . To 
make a well-defined situation, we assume a finite system with periodic boundary conditions 
having dimensions L x x L y where L x is even and L y is odd. This forces an odd number of 
domain walls (which must be stripes) running in the x direction. These domain walls are 
composed of 1/2 of a hole per column and we call them stripes. This model has been studied 
by ED for spinless fermions and hard-core bosons, and the stripes were shown to be stable 
against phase separation, for fermions When the number of particles is L x (L y — 1)/2, 

only one stripe is allowed. If we remove a few more particles (less than L x ), it is energetically 
favorable for the holes to attach to the stripe, since holes off the stripe can only form confined 
droplets [5]. From here on, we reserve the term "hole" for additional holes beyond those 
needed to create a stripe. In the case of a single undoped stripe, the boson and fermion 
models possess the same energy spectrum [2]. 

On an undoped stripe, hops of particles are equivalent to stripe fluctuations (see fig. Q. 
If we define the stripe height y(x) to be the mean of the y coordinates of the first particles 
above and below the stripe, in column x, we can map the up and down steps of the stripe 
height to "spins" in one dimension by defining s(x) = \[y{x) — y(x — 1)], taking values ±1/2. 
The corresponding Hamiltonian is the one-dimensional spin-1/2 XY Hamiltonian 

We, = where H\ x = S+ S m + <?r S+ +1 , (2) 

i 

where S 1 - , S~ are spin raising and lowering operators. This model can be solved by the 
Jordan- Wigner transformation, in which we replace each up spin by a spinless fermion, and 
each down spin by an empty space. The Hamiltonian is replaced by a non-interacting hop- 
ping Hamiltonian for the fermions. Thus, the ground state of a horizontal, undoped stripe 
is equivalent to that of a (one-dimensional) half-filled sea of free (spinless) fermions with dis- 
persion e(fc) = — 2cos/c. From this, the ground state energy, for large L x , is 6q L ^ = —2L x /-k. 
This implies the chemical potential of the stripe is fj- s tripe = — 4/-7T, since half of a particle 
is removed per unit length of a stripe. In a more general case, we may force a stripe with 
some overall tilt by replacing our rectangular boundary conditions with (L x , b) x (0, L y ), for 
L x + b even. In that case the spin chain has total spin 6, mapping to ID fermions results in 
N =(L x + b)/2 fermions and the ground state energy is e 1 ^*'^ — —2L x sm(Nir/L)/ir |2"lll(J|. 
Note that a diagonal stripe (b — L x ) has zero energy, as there are no allowed hops. 

One hole on a stripe. - If an additional hole is added to an a stripe, the energy is lowered 
due to the hopping energy of the hole along the stripe. We define the ground state energy 
difference A = E^oie — E str i P e, where Ehoie, E str i pe = are ground state energies of a 

stripe with a hole and an undoped stripe, respectively. For stability of a stripe array state, 
we must have A > Stripe ; otherwise, doping would add holes to existing stripes, probably 
forming a phase-separated droplet. To obtain a variational bound for A, we will make (below) 
an approximation by using a subset of the Hilbert space. 
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Once a single hole is added to column x of the stripe, the stripe height difference in the two 
columns adjacent to the hole can take values of or ±2 rather than ±1, i.e. s(x), s(x+l) 6 
{0, ±1} (see, for example, figs. |21 • Thus, we get a pair of "spin-1" impurities on both sides 
of the hole, at positions x and x+l on the spin chain. This looks like a Kondo model with 
two mobile impurities, that are bound together. Similar systems, with t-J Hamiltonian and 
mobile spin-1/2 impurities in a spin-1 chain, have been the subject of theoretical study, as a 
model for the charge transfer insulator Y2_ a: Ca a ;BaNi05 11 . However, our model is quite 
different from those, because the hopping matrix is quite elaborate, reflecting the allowed 
moves in the original stripe. 

In order for our variational scheme to work, we need to separate the "hole degrees of 
freedom" from the "stripe degrees of freedom". To make this separation apparent, we find 
it more convenient to introduce a different, equivalent notation, in which the stripe and the 
hole is represented by a one dimensional spin-1/2 chain of length L' x = L x — 2, with three 
additional "particles": a hole (represented by a dot), that marks the position of the hole in 
between two spins, and two "brackets" (right and left), each taking the place of a spin. If 
s(x+l) = 1 (—1), then the right bracket takes the places of the first down (up) spin to the 
right of the hole, and s(x+l) is set to 1/2 (—1/2). If s(x+l) = 0, the right bracket is placed 
in position x+l. Similarly for the left bracket and s(x). Note that the total spin is preserved 
by this mapping. For example, if we use a double arrow to denote s(x) = ±1, 

ITIII-U-'frllt => I LI I I I • TJ I T 
IT44|OJ||tt => I T I I IU IJ T (3) 
tlttTftltttl =*• T LT T T T • T T TJ 

If stripe fluctuations occur in the vicinity of the hole, it might become stranded, i.e. isolated 
from the stripe, as in fig. El However, we observe that in the ED ground state of ref. 0, the 
probability for this to occur is negligible (about 0.03 for an untilted stripe and less than 0.01 
for a tilted one). Thus, in the following discussion we suppress hops that lead to the stranded 
state. Under this assumption, the spin chain is free to fluctuate only outside the brackets. 
Each of the three additional "particles" can hop by one step freely, as long as all of the spins 
between the hole and each of the brackets are in the same direction, and as long as their order 
is preserved (i.e. they do not hop across each other). This mapping can be shown to be exact. 

The Hamiltonian for the hole on the stripe in this model can now be broken up TLhoic = 
TLh + H r + TLi + H' sp , where Hh, H r , Hi are the hopping terms for the hole, right bracket 
and left bracket, respectively and H' sp is the spin exchange Hamiltonian acting only outside 
the brackets. The Hamiltonian TChoie acts on the Hilbert space of "allowed states", i.e., all 
states in which the spins between each of the brackets and the hole are in the same direction. 
We denote the state of the hole and two brackets by \x,l,r), where the hole is between the 
spins at positions x and x+l, there are / spins (sites [x — l +l,x]) between the hole and left 
bracket, and r spins between the hole and the right bracket (sites [x+l, x+r]). For example, 
in the three examples in above, if we number the spins starting from 1 at the left, the 
hole+bracket states are 5,4,1), 5,0,2), 5,4, 3), respectively. Using this notation 

H h = - ^ ^>o (^f 5 + Ml ~ 5 )) \x,l,r)(x + lj + l,r-l\+h.c, (4a) 

H r = -Y,( V ^r+^rA^~'Px+r))\xJ,r+l)(x,l,r\+h.C., (4b) 

Hi = -Y J {'Px-i + W- r P^ l ))\x,l + l,r){x,l,r\+h.c, (4c) 
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Fig. 1 - Particle hops (stripe fluctuations) in an undoped stripe are equivalent to spin exchange. 

Fig. 2 - Hole stranding: Stripe fluctuation across the hole can leave the hole away from the stripe, 
with no possible allowed hops. We ignore these moves our variational model. 

Fig. 3 - Three successive moves that shift the hole by one position to the right. The sequence shown 
here is a main path sequence that can occur regardless of the configuration of the stripe. 



where V^ s = 2S , |5| +1 + 1/2 is 1 if S* = S* +1 and zero otherwise, and denotes a sum 
only on sites that are outside of the brackets (i ^ [x — l,x+r]). We have thus obtained that, 
using the mapping, the coupling between the hole and spin degrees of freedom seems to be 
limited to the interval between the brackets. Therefore, it is natural to construct a variational 
wavefunction by starting with a complete decoupling between the hole and the stripe, and 
project out the "illegal" states, in the spirit of the Gutzwiller projection. 

Motion of a hole on a stripe. - In order to examine the motion of a hole on a stripe 
we find it instructive to represent the various states of the stripe by a state graph [5] , which 
is defined such that each node represents a state and each line represents an allowed hop 
between states. In order for the hole and brackets to shift by one position, in this model, 
three individual hops are required (see, e.g., fig. [3J|. Generally, it is clear that the allowed 
hole and bracket moves depend on the background stripe configuration. The simplest possible 
motion is when each of the brackets does not stray more than one position away from the 
hole, thereby eliminating any dependence on the stripe. This motion can be represented in the 
state space by the solid lines in fig. 01 and we refer to it as the "main path" . Depending on the 
local tilt in the stripe (runs of spins in the same direction), the available states are increased. 
In the extreme case, if the stripe is a 45° diagonal, then the hole and brackets can be treated 
as three non-interacting fermions and their energy is lower than in the untilted case. This 
implies that the hole prefers the stripe to be locally tilted around it. However, the spin chain 
loses energy by having a region where all of the spins are the same, since no spin exchange is 
possible there. We would expect the interplay between these two competing effects to result 
in a small tilted region ("kink") around the hole. This has indeed been inferred from ED of 
ref. 0, and is one of the phenomena that we look for in our solution for this model. 

Examining the motion of a single hole on a stripe, we find that any sequence of moves that 
returns the system to its original configuration involves an even permutation, except for hole 
hopping around odd boundary conditions This implies that for even or infinite L x , the 
boson and fermion spectra are identical. 




(4d) 
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Fig. 4 - State space representation of hole motion for a particular background stripe/spin configu- 
ration. The 1 symbol represents an unspecified spin (j or j). The solid lines represent the "main 
path", in which the hole motion is decoupled from its environment. The dashed and dotted lines 
represent the additional available hops when there is a run of two or three spins in the same direction, 
respectively. If we read this from left to right, a horizontal line represents one hop of the hole to the 
right, a downward(upward) diagonal line represents hopping of the right(left) bracket. 



Variational wavefunction. In order to calculate the hole spectrum variationally, we 
introduce a projection operator V^f{x) acting on states of a spin chain of length L' x , to the 
sub-space where all of the I spins at sites [x— x] are in direction a, and r spins at positions 
[x+l,x+r] are in direction (3 (a, (3 G {!,?})■ T^if^) annihilates states that are not in this 
subspace. We define an orthonormal basis set: 

xf(x) = j^\x,l,r}V? r f} (x)\* F ), (5) 

Ir 



where is ground state of an unrestricted spin- 1/2 ring of length L' x ; M^f = yj (V"f) F is 
a normalization factor, where (V^P)f = 'f\V^ f) is independent of x. 

Now we want to calculate (xff '(x)| Whole Ixl-v ( x ))- The matrix elements for Hh, H r , Hi 
are straightforward, but those of H' sp are a little harder to calculate. Since H' 8p does not 
change the hole+brackets state |x, I, r), we only need matrix elements between pairs of states 
with the same x,l,r, i.e. (xfj 3 ' (x)\'H' sp \xf r ^ ( x ))- These can be calculated using the fact that 
Tt' commutes with Vy{x) 

(xf (x)\K\xt7(x)) = ^y ^wE'^^w)^ 



\'lr IF 



(6) 



which is independent of x. We used: (i) \H\ X , V"f(x)] — for i £ [x— l,x + r] (modulo L' x ). 
(ii) P"f(x)W ex — for i e (x— l,x) or i S (x,x+r). Note that when 2 = 0, a is undefined, so in 
order to generalize the notations, we arbitrarily set a =J, when 1 = (Vjf r = 0), and similarly 
for r,(3. We defined 8 l = S a ^ if I > and 5' = otherwise. In order to calculate the 
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Fig. 5 - Circles: Ground state energy difference A = Eh i e ~E str ipe as a function of overall stripe tilt, 
for infinite L x . The dashed line represents the chemical potential jj, st ri P e- {fJ-strive < A is required for 
stripe stability). Squares: Main path probability P m am- 

Fig. 6 - Energy difference A obtained in ED and variational calculations, for finite system sizes. 
Results are shown for horizontal stripes (b = 0), as well as tilted boundary conditions (b = 1,2). 
The solid lines are linear interpolations of the ED data. Note that although A decreases as the tilt 
increases, Ehoie generally increases with tilt. 



energy spectrum using the variational wavefunction, we Fourier transform xtr ( x ) an d nn d 
{xtr ' (k)\T~Lhole\xTr^ s )) m a straightforward manner. 

Results. - The results presented here were obtained for a subset of the variational basis 
set with l + r < 10 for our calculations. If one does not take advantage of mirror symmetries 
between states, the Hamiltonian is a 221x221 matrix. We find the ground state and ground 
state energy of this Hamiltonian as a function of L x , b, and k. Increasing the variational basis 
set to up to l + r < 13 results in a relative change of less than 0.1% in A for the maximum 
overall tilt presented here, and less than 10~ 14 for zero tilt. In our calculation, we take into 
account the fermion statistics only by adding a phase to hole hops in (finite) odd boundary 
conditions. This is justified because these are the only moves that can induce an odd cyclic 
permutation of the particles ^2] ■ 

Figure [S] shows the ground state energy difference A, as a function of the overall stripe tilt. 
We obtain A = —0.29 for an (infinite) untilted stripe. As the stripe tilt is increased, both 
the energy of the undoped stripe, E str i P e, and Ehoie increase, because stripe fluctuations are 
reduced. However, the difference A decreases, because the hole's kinetic energy is enhanced 
by the additional tilt. Figure HO shows A for finite system sizes, compared to respective ED 
results of ref. [2], and their extrapolation to L x — > oo. The variational energy is, of course, 
an upper bound to the true energy of the system. In the ED of ref. 0, Ehoie was minimal 
for stripes with slightly tilted boundary conditions (b — 2 for even L x , 6 = 3 for odd L x ). 
This suggested that the hole tends to bind to a kink in the stripe in order to increase its 
kinetic energy, at the expense of stripe fluctuations, forming a polaron. In our variational 
calculation, we observed the same effect for small system sizes (up to L x KilQi), however, for 
larger systems, the minimum energy is for b = 0(1) for cven(odd) L x . This indicates that 
the preference for tilted boundary conditions, as observed in ED, may be only a finite size 
effect. We also calculated the hole dispersion. On an untilted stripe, we find an effective mass 
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Table I - Variational and ED ground state probabilities for (10, 0) x (0, 7). All other states, except 
those related by symmetry to the ones presented have probability of less than 1 %. 



Basis State 


Xoo 


it 
Xol 


Xn 


Xn 


X02 


variational 


23.4% 


13.7% 


4.3% 


2.9% 


1.4% 


ED 


19.2% 


12.3% 


4.7% 


3.6% 


1.8% 



m* h = 6.43m*, where m* — Q5 is the mass of a non-interacting particle hopping on the lattice. 
This agrees remarkably well with ED, which gives m* l ~Q.7m* |2J. If a hole would be forced 
to move only on the "main path", its effective mass would be 7.35m*. When the stripe has 
an overall tilt, the hole mass is reduced significantly, e.g. at dy/dx = 0.5, m^ — 4.57m* . 

A good measure of the coupling between the hole and the stripe is p ma im the ground state 
probability of being in one of the "main path" basis states. In the variational ground state 
for an untilted stripe, we find p m ain = 0.91, i.e. the hole tends to be decoupled from the stripe 
configuration. This probability is reduced as the stripe tilt increases and it is easier for the 
brackets to move from the hole (see fig.|SJ). This trend is supported by ED calculations, but 
the variational values for p ma i n are higher by up to 0.1. Comparison of the variational and ED 
ground states, reveals that the our calculation does a good job of qualitatively capturing the 
composition of the exact ground state from the basis set (e.g. table |]J, but it overestimates 
the weight of states with low l+r. 

In summary, we calculated the energy and the ground state for a single hole on a stripe, 
using a mapping to one dimension and a variational wavefunction constructed to decouple the 
hole and stripe degrees of freedom. We did not find evidence that the hole binds to a kink in 
the stripe, in the untilted case, for L^^IO. 

* * * 
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